General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



UNIV£4> 


N83-22<i97 


(NAS&-CB-170163) IHELLIS PHASE CODES EOE 
POHEB”BAMDHITH EFFICIENT SATELLITE 
COHBUHICATIOHS Fiaai Report, 10 Apr. 1979 - 
22 Feb. 1981 (Virginia Univ.) 79 p Onclas 

HC A05/HF A01 CSCL 17B G3/32 15260 


Final Report 
NAS5-25634 

TRELLIS PHASE CODES FOR POWER/BANDWIDTH 
EFFICIENT SATELLITE COIVIMUNICATIONS 

Submitted to: 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland , Ohio 44135 


Submitted by: 

S. G. Wilson 
Principal Investigator 

J. H. Highfill 
Senior Scientist 

C-D Hsu 

Research Assistant 

R. Harkness 
Research Assistant 





COMMUNICATIONS SYSTEMS LABORATORY 

DEeiARTMENT OF ELECTRICAL ENGINEERING 
SCHOOL OF ENGINEERING AND APPLIED SCIENCES 
UNIVERSITY OF VIRGINIA 


& 


Final Report 
NAS5-25634 

TRELLIS PHASE CODES FOR POWER/BANDWIDTH 
EFFICIENT SATELLITE COMMUNICATIONS 

Submitted to: 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland , Ohio 44135 

Submitted by: 

S. G. Wilson 
Principal Investigator 

J. H. Highfill 
Senior Scientist 

C-D Hsu 

Research Assistant 

R. Harkness 
Research Assistant 



Department of Electrical Engineering 
RESEARCH LABORATORIES FOR THE ENGINEERING SCIENCES 
SCHOOL OF ENGINEERING AND APPLIED SCIENCE 
UNIVERSITY OF VIRGINIA 
CHARLOTTESVILLE, VIRGINIA 


Copy No . 



/> 


Report No. UVA/525634/EE81/101 
April 1981 


1, Report No, 

2, Government Accession No, 

3, Recipient's Catalog No. 

4. Title and Subtitle 

Trellis Phase Codes for Power/Bandwidth 
Efficient Satellite Communications 

5. Report Date 

Anrll 1981 

6, Performing Organization Code 

NAS 5-25634 

7, Auihor($) 

S. G. Wilson, J. H. Highfill 
C-D. Hsu, and R. Harkness 

8. Performing Organization Report No. 

UVA/525634/EE81/101 

10. Work Unit No. 

9. Performing Organization Name and Address 

Communications Systems Laboratory 
School of Engineering and Applied Science 
University of Virginia, Thornton Hall 
Charlottesville, VA 22901 

11. Contr ct or Grant No. 

I^.Tyi^ ^ Report^and Period Covered 

04/10/79 - 02/22/81 

12. Sponsoring Agency Name and Address 

NASA Lewis Research Center 
21000 Brookpark Road 
Cleveland, OH 44135 

14. Sponsoring Agency Code 


15, Supplementary Notes 


16. Abstract 


This report summarizes work performed in support of improved power and spec- 
trum utilization on digital satellite channels. Specific attention is given to the class 
of signalling schemes known as continuous phase modulation (CPM). 

The specific work described in this report addresses 

(1) analytical bounds on error probability for multi-h phase codes, 

(2) power and bandwidth characterization of 4-ary multi-h codes, and 

(3) initial results of channel simulation to assess the impact of band-limiting 
filters and nonlinear amplifiers on CPM performance. 


17. Key Words (Suggested by Author(s)) 

modulation and coding 
satellite communication 

18. Distribution Statement 

19 Security Oassif. (of this report! 

Unclassified 

20, Security CUssif (of this page! 

Unclassified 

21, No. of Pfiges 

79 

22. Price 




Foi sale by the National Technical Intotriation Seivice, Springfield Virginia 22161 



TABLE OF CONTENTS 


1.0 INTRODUCTION 

2.0 ERROR BOUNDS FOR MULTl-H PHASE CODES . . . . . 

I. INTRODUCTION 

II. BACKGROUND 

III. DEVELOPMENT OF ERROR BOUNDS, UNLIMITED 

PATH MEMORY 

IV. UPPER BOUND FOR FINITE MEMORY 

V. NUMERICAL EXAMPLES . . . . 

REFERENCES ..... . 

3.0 4-ARY MULTI-H CODES 

REFERENCES ........ 

4.0 CHANNEL SIMULATION 


4.1 SIMULATION RESULTS . . , . 

4.2 SOFTWARE DESCRIPTION . . 


REFERENCES 


APPENDIX A 


1.0 INTP.ODUCTION 

In this final report on "Trellis Phase Codes for Power/Bandwidth 
Efficient Satellite Communication," (NAS5-25634) we summarize work 
performed since the last report of July 1980. Since that time we have 
continued our study of multi-h and partial response FM digital modula- 
tions, and have emphasized three special aspects: 

1. analytical error bounds for multi-h phase codes 

2. a study of 4-ary multi-h codes 

3. channel simulation for partial-response and multi-h codes to 
assess the effects of handlimiting and nonlinear amplification. 

These are discussed in detail in following chapters of the report. 

Additional work performed under the contract has been described 
in two pervious semi-annual reports to NASA/Lewis Research Center, 
published in November 1979 and July 1980. 
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2.0 ERROR BOUNDS FOR MULTI-H PHASE CODES 


This section is a revision of a manuscript accepted for publication 
in the IEEE Transaction on Information Theory. Numbering of equa- 
tions and references is internally consistent. 

I. INTRODUCTION 

Multi-h phase codes, described in detail by Anderson and Taylor 
[1], and earlier conceived by Miyakawa, et al. [2], represent a class of 
constant envelope signal designs providing attractive gains in the 
power/bandwidth tradeoff, relative to MSK or QPSK. Time-variation of 
the transmitter frequency deviation parameter among a small set of 
rational numbers provides delayed remergers in the phase trellis. This 
in turn provides increased minimum distance between signals which may 
be exploited by a maximum likelihood sequence estimator to enhance 
detection efficiency. Analysis thus far of these signals has concentra- 
ted on minimum distance to yield asymptotic behavior, and some simula- 
tion results are available, [3,4]. 

These signals may be represented in terms of a finite-state trellis, 
though a time-varying one due to the cyclical nature of the modulator 
deviation parameter. In order to determine bounds on error probabili- 
ty, we extend the union bound/transfer function approach earlier used 
for linear convolutional codes [5] to the multi-h phase code case. In 
closely-related work, Aulin [6] has developed similar bounds for the 
class of partial-response FM codes with fixed modulator deviation. As 
with earlier applications, the principal difficulty is in formulating the 
state diagram and the necessary branch gains. Once this is accom- 
plished, numerical evaluation is straightforward. 
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Lower bounds on receiver performance are more easily obtained in 
terms of a single (or few) neighboring paths, and it will be seen that 
the two bounds are close for performance levels of typical interest. 

The paper is organized as follows. Section II provides a brief 
review of the multi-h signal structure, and its state-space description. 
Section III develops the upper and lower bounds for the case of unlimi- 
ted path memory, incorporating the concept of difference-state trellises 
and the associated transfer functions. The effects of finite decoder 
path memory are addressed in Section IV. The paper concludes with 
several illustrative examples comparing the analytical bounds with simu- 
lation data . 

II. BACKGROUND 

We assume the transmitted signal is of the constant-envelope form 

s(t) = (2?)"^ cos (Wj,t + Qq + <l>(t,a)) (l) 

where P is the power of the symbol, 0^ is an arbitrary phase angle, 
and (|)(t,a) is the (causal) phase modulation induced by a semi-infinite 
sequence 


3 j ) ^2 * ■ * ‘ ^ 

The symbols a^^ are selected independently each T seconds with equi- 

probability from the set {-(M-1), , -3, -1, 1, 3, .... (M-1)} so 

that the signalling is M-ary. 

Over the n^^ signalling interval, (n-1) T < t < nT, the phase term 
in multi-h coding is a continuous waveform given by 
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(|)(t,I) = n I a. h. + n h g(t - (n-l)T)dt (2) 

j=l J J (n-l)T “ 

The first term in (2) gives the accumulated phase due to symbols 
through while the second term provides the time-varying phase 

increment over the n^^ interval. Also in (2) we have defined g(t) to 
be a frequency pulse for spectrum-shaping purposes. We assume that 
g(x) is non-zero for 0 < t < T, and that 

;Jg(t)dt = 1 (3) 

Under this normalization h^^ plays the role of a modulation index, 
and hj^/T is proportional to the average frequency deviation over the 
symbol interval. We remark this description is more general than that 
of [1], in that M-ary modulation, as well as nonlinear phase trajecto- 
ries, are allowed. 

In standard digital FM transmission, h^^ is some constant h. For 
MSK designs, h = 1/2, and for binary CPFSK, h = 0.715 maximizes the 
detection efficiency. However, in multi-h coding, h^^ is cyclically 
chosen from a set of K rational numbers: 


Pi P2 


«K= ^ T’ f’ ^ ^ 


q 


.(1) u(2) 


(4) 


That is 


and thereafter 


hj = hj^ = 


h . = h 

n+jK n 

The set of possible phase trajectories associated with data sequen- 
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ces of length n is a phase tree of trajectories. By virtue of the 
rational“index assumption, it will be seen that the number of distinct 
phases, modulo - 2n, attained at the end of symbol intervals is exactly 
q. (A complication is that there may be a different q values for differ™ 
ent times, but this does not influence our development). Thus, the 
tree collapses to a phase trellis, describing the allowable phase transi- 
tions, and tlie allowable signal trajectories. 

By defining the signal state at time n as 

n-l 

S - n Z h. a., mod lit (5) 

n j=i J J 

we obtain a fihite-state description of the signal, whereby the state 
evolves among q states (phases) according to 

S , , = S + n h a , mod 2n (6) 

n+1 n n n’ ^ ' 

This description is essential for maximum likelihood decoding via the 
Viterbi algorithm. 

As a unifying example throughout, we shall consider a simple 
binary code defined by Hg = {2/4, 3/4). The state trellis for this 
example is shown in Figure 1, where it may be seen that the allowable 
state transitions are time-varying, depending on which deviation is 
currently in effect. The trellis shows state transitions, and is not 
meant to indicate actual phase variation, although in the usual case the 
phase segments are piece-wise linear. The trellis is structurally very 
regular, its appearance complicated only by the mod 2n phase defini- 
tion. 

For this example, it may be seen that two data sequences which 
differ in the first position produce state sequences which are unmerged 
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for at least K + 1 = 3 bits in the trellis. Generally speaking long 
remerger lengths are desirable for energy efficiency, and it may be 
shown that there exist sets Hj^ which achieve an unmerged span of K+1 
levels provided q > M , [1), [4]. 

Assuming coherent reception of the multi-h signal in white 
Gaussian noise, the optimal sequence detector is comprised of a bank of 
correlators coupled with a q-state Viterbi algorithm decoder. The 
correlators in effect generate branch metrics for the trellis search. 
The receiver locates the data sequence producing the largest waveform 
correlation, ignoring for the moment questions of quantization and path 
memory. An alternate view is that the receiver does minimum distance 
decoding, with distance measured in the L 2 norm. Distance between 
waveforms is key concept in the development, and is given by 


■=12' = -Ir ^ 

N 


The '.lormalization is selected so that the probability of error in a two- 
signal test is 


P 

e 


= Q 



1/2 


( 8 ) 


' 0 ' 

where E = PT is the energy per symbol, and is the one-sided noise 
spectral density. With this definition, antipodal signals of length T 
have a normalized distance of 2. 

Quantities of special importance are the minimum distance over an 
interval of length N, defined as 


d . 


nun 


2 

N 


min min 

i j 



C9) 
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and its limiting value in N, d£j,gg. These will determine the detection 
efficiency at large E/N^^. 

III. DEVELOPMENT OF ERROR BOUNDS, UNLIMITED PATH MEMORY 
In this section we proceed to develop upper and lower bounds on 
the probability that a decoding error event begins at some level n, 
denoted P^, and on the probability of symbol error, denoted P_. The 
receiver memory is assumed arbitrarily large. The development utilizes 
standard transfer function bounding, e.g., [5], modified to incorporate 
difference state concepts [6] and the time- varying nature of the trellis. 

An error event uegins at time n if the decoder selects sequence 6 
rather than the transmitted sequence a , where a and E first differ 
at time n. The probability that two sequences are confused depends 
only on their phase difference trajectory (as well as E/N^). Since 
phase is in turn proportional to the data sequence, we need only consi- 
der possible data difference sequences, rather than consider all possible 
sequence pairs (a , E ). This fact has been utilized by Aulin [6] for 
the case of partial-response FM trellis codes. 

For multi-h codes, the difference-state for a pair of sequences is 
given by the accumulated phase difference, mod 2n. That is 

n-l 

A = I n h. Y- > mod-27t (10) 

^ j=l J 

where Yj = Sj " hj. The difference-state trellis shows the possible 
evolution of difference-state among the q possible values for (10). 
Since Yj is in a set of size 2M-1, there are 2M-1 branches entering and 
leaving each node of the difference-state trellis. Figure 2 illustrates 
the difference-state trellis for the {2/4, 3/4} binary code described 
earlier . 
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Error free decoding corresponds to residence in the * 0 state, 

i.e. the selected path agrees with the transmitted path. An error 

event begins at time n if the selected path splits from the zero differ- 

til 

ence- state at the n level. The probability of this event, P , is 
upper bounded by the sum of probabilities that the decoder selects a 
specific difference sequence which leaves the zero difference-state at 
time n and remerges later in the trellis. 

Actually this bound is unduly pessimistic since transmitted sequen- 
ces are not in one-to-one correspondence with difference sequences. 
For example consider M = 2. The allowed values for Vj^ 2, 0, -2. 

= 2 can be' formed only by (a^^, b^^) = (1,-1), whereas can 

correspond to (a^, b^^) = (1,1) or (-1,-1). To properly average with 
respect to transmitted sequences, we must give lesser weight to the 
influence of = ± 2 transitions. Similar effects happen with larger M, 
and Aulin [6] has shown a multiplying factor 

c = n [M - |y^|/2) /M (11) 

n 

must be applied to the probability of each difference sequence to pro- 
vide the proper averaging. 

Assuming for the present that deviation h^^^ is effect at time n, 
we can now write the conditional event error probability as 

» -d.^ E 1/2 

'ell" 

1=1 0 

where i is an index on the set of unmerged difference sequences, jd.^ is 
the distance corresponding to the i*^" sequence conditioned by the fact 
that and c. is a weighting factor obtained via (11). 
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Signal flow graphs and their transfer functions may be used to 
evaluate (12) numerically. The difference-state transition diagram is 
viewed as a flow graph with zero state split into input and output 
nodes. Paths are labelled with a gain of the form cD^I^L where c is 
given by [M - |yI/2]/M, x is the distance accumulated in undergoing a 
given transition, and y is an error indicator (y = 1 if y 5^ 0 and y = 0 
if y =: 0). 

In the multi-h case, the state transition diagram is time- varying. 
Whereas a single transition matrix suffices in the usual case, K state 
transition diagrams (and matrices) are needed here. Returning to our 
example, we show the two split difference-state diagrams for the {2/4, 
3/4} code in Figure 3. In tracing paths, one must alternate between 
diagrams. 

To handle the time-varying case, we imagine an ensemble of trans- 
mitter/receiver pairs for which the deviation index at time n = 1 is 
randomly selected among K values equiprobably. Thereafter each 
system runs according to the usual mulLi-h format. One then sees that 
the error probability at any time n is an average of error probabilities 
conditioned on a specific choice of deviations, 


P 

e 




(13) 


and the same holds for symbol probabilities. More practically, one may 
view this in terms of a time average for a single realization. 

Next, we upper bound the error probability noting the other 

Pgjj follow directly. We define, Tj^ to be the matrix of path gains 
between difference states when deviation h^^'^ is in effect, T 2 to be the 
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matrix corresponding to and so on. Each matrix will be square 

with dimension q + 1, arising from q difference states plus one for 
splitting of the zero-d’fference state. We index the difference states 1, 
2 , . . . , q , and let 1 and q + 1 be the indices for the split zero state . 

-k,T. 

The jk*"” element of represents the path gain from state j to state k, 
and will be a polynomial in D, I and L. Figure 3 itemizes these matri- 
ces for the example code. 

In finding the transfer function of the graph, we note that 
represents the possible actions after a single step, T^^T 2 represents the 
actions after two steps, etc. After any prescribed number of steps, 
the upper right-hand corner elemefit of the sum matrix 

4 Tj + T^T2 + T^T2T3 + . . . . (14) 

provides the path enumeration we seek, namely the set of possible path 
distances, their lengths, and respective number of symbol errors. 

In the case of unlimited memory, we have an infinite sum of matrix 
products in (14). Noting the periodic nature of the h^^, (14) may be 
regrouped as 
00 

= [ Z (TjT 2 ... Tj^)-^] [Tj + TjT 2 + ... .. Tj^] 

j=0 

= [I - T^T2 ... [Tj + T^T2 + ... T^T2 ... (15) 

We call the upper-right corner element of the transfer function of 
the graph, when h*'^^ is the initial deviation, and denote it Gj^(D,I,L). 
It may be written in general as 

00 2 . p. 

G (D,I,L) = Z c. D^i (16) 

^ i=l ^ 
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where Cj is a positive constant, is a distance associated with an 

error event of length A., v. is the corresponding number of symbol 

errors. There may be many error events of a given length having 
2 

differing d. and v^. 

Recalling the upper bound to event error probability of (12) we 
see that if the series expansion (16) could be obtained, the bound is 
obtained. Such is not possible except for simple examples, so we seek 
an efficient numerical procedure. 

By the inequality [5] 

Q (Ca+b)^'^^) < a, b > 0 (17) 

we may write (12) as 

2 2 
„ I/O ■ E » d.'^E 

En1/2n rl mxn. ^ \ /*io\ 

— ) ) exp [— ~2 n ^ ^ c^expC- — ) (18) 

o i=l o 

But the summation in (18) is precisely Gj^(D,I,L) with D = e’^'^^^o, I = 
1, L = 1. Thus, conditioned on h^^^ being n^^ deviation, we have 


el 


< Q(0 


id^. 

1 mia 
N 


2 2 
,d^. E -d'^. E 

Tk I >> r\ /* / ^ mxn s i./2\ / j. rnxn \ ^ ^ 

^e 1 ^ N ^ ^ ^~2N ^ G^(D,I,L) 

o 0 


_ -E/2N 

D = e o 

I = L = 1 

(19) 


< G^(D,I,L) 


2n ,d . E , 

T» /o\T / / 1 iTixn ^ x/2 

D = f ii > 

I = L = 1 / ° 


The latter follows by the inequality 

Q(x) < e"^ ^^/(2n)^'^^x 

which weakens the bound slightly for small signal-to-noise ratio. 


(20) 


15 


Q. 


POOR quality 


Furthermore, the conditional upper bound on symbol error probability is 



tVi 

where is the number of symbol errors associated with the i ^ error 
event. Using steps as above, this becomes 


P 

s 



3Gj(D,I,L) 

91 




The desired derivative may be numerically evaluated using finite differ- 
ences . 

To complete the analysis we simply must average these bounds with 
respect to the K possible values for h^^\ as in (13), and the final 
upper bounds become 


1 ^ 
P < 4 I 
e K , . 
1=1 


2rc 


.d^. E 
1 min 

N 


■ 1/2 


G^(D,I,L) 


-E/2N 

D = e ' o 
I = L = 1 


(23a) 


- 1/2 


P < ^ 2 2n 
s K . 

1 


,d . E 
1 min 

N 


9G.(D,I,L) 

ai 


n -E/2N 
D = e ' o 

I = L = 1 


(23b) 


Lower Bounds 

Although upper bounds on performance are generally of most 
interest, and in fact it is known that the upper bounds derived here 
are asymptotically tight for large E/N^, lower bounds may be of interest 
and are easily calculated. One principal use is in assessing the region 
for which the upper bound is tight. 
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A lower bound on error event probability is provided by computing 
the probability that any error sequence diverging at level n has higher 
metric (correlation) than the transmitted sequence, and we can make 
the lower bound tightest by choosing that neighbor sequence having the 
minimum distance. As before, there are K equiprobable selections for 
the n " deviation, and a simple averaging of conditional lower bounds 
provides the desired lower bound. The minimum distance depends on 
the initial deviation. This minimum distance may generally be found by 
listing short remerger events of length L defined by 

L 

I h Y =0, mod 2n (24) 

n n ' 

U”1 

and' evaluating the distance from the difference state transition dia- 
grams, In complicated cases, d^.^^ may be found using dynamic pro- 
gramming on the difference-state trellis. (This is the same 
needed for the upper bound). 

These lower bounds to error event and symbol error probabilities 
become 

1 K .d^. E , 

> i I c. Q ) 1/2) (25,) 

i=l “o 

2 

^ I ^ ''i '"i Q(( - - |— ^ ^^^) (25b) 

i=l o 

where v. is the number of symbol errors on the various minimum dis- 
tance paths. The bounds of (25) use only the minimum distance event 
for each possible starting deviation, and one finds the bound is weak at 
low signal- to-noise ratios, particularly so when several error events 
have comparable distances. 
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IV. UPPER BOUND FOR FINITE MEMORY 

The infinite path memory assumption of the previous section is of 
course not realizable in practice, although one expects that if the 
decoder delay is suitably large, then ideal performance is closely ap- 
proached. Since we would like to operate with minimal delay and memo- 
ry, questions arise about the effect of a given delay on performance. 

We assume the decoder operates as follows. The Viterbi algorithm 
performs trellis computations as usual, but maintains survivor paths in 
memory for only N symbols. Thus at level n + N of the trellis, the 
path with the currently best metric is located, and the oldest symbol on 
its path map, corresponding to a^^ is released by the decoder. Trellis 
pruning is not performed if survivor paths disagree in the oldest posi- 
tion. Relative to the case of pruning of such paths, P_ actually is 
smaller, and the analysis is simplified. 

We classify error events into two classes: those made by an ideal 

ML decoder operating with infinite memory; and those due to truncation 
wherein an unmerged path has highest metric at the current time, and 
the oldest symbol in its history is released by the decoder. Following 
Hemmati and Costello [7] , we sum the error probabilities for each type 
as an upper bound on overall error probability. This allows some 
double-counting of error events. 

Consider again the difference-state trellis, where the all-zero 
sequence corresponds to error-free operation. An error event due to 
truncation can occur at time n whenever a path splits from the zero 
difference- state at or before time n and remains unmerged at time n + 
N. With the aid of Figure 4, we illustrate error events of both types. 
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Events 1 and 2 are of the type which could produce symbol errors at 
time n in ML decoding, and were included in the analysis of Section III. 
Event 3 diverges at time n and remains unmerged at time n + N, so it 
leads to a truncation error at time n. Finally, path 4 is a sequence 
which diverged earlier, remains unmerged at n + N, and may produce a 
symbol error at time n. 

We assume for notational convenience that the decoder truncation 
depth N is a multiple of the cycle length K, i.e. N = rK. Since multi-h 
codes have apparently little theoretical or practical benefit for K > 4, 
this represents no real instrumentation constraint. 

As before ’ we analyze performance for a specific deviation condi- 
tion, then average with respect to the K such choices. We first let 
( 1 ) 

h^ ' be in effect at the time of interest, n. We wish to enumerate all 
difference sequences which split from the zero difference-state at time n 
or earlier, and which remain unmerged at time n +• N. Consider the 
matrix (Tj^Tg.. This enumerates paths of length N = rK, and 

the elements in the top low, except the upper-right corner element, 
constitute the unmerged paths which split from the zero state at time n. 
Similarly, Tjr(T^T 2 -. T^r)^ describes paths of length N + 1 diverging at 
time n - 1. Thus we sum matrix products of order N or larger, and 
define this sum to be 1^^: 



[ 

[ I + 


" Tk-iTk 
* Tr-iTk 


+ . . . 


* TjTj .. Tg + 

T2T3..T^1 [ I - 

(26) 
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Finally we let (aj^(D,I,L) be the polynomial sum of all top-row entries in 
excepting the corner element, and term this the transfer function 
for truncation error. Following the development of Section III, we have 
that, conditioned upon h^^^ being in effect at time n, the symbol error 
probability due to truncation is upper-bounded as 


S 


trunc 1 






1/2 


n - -E/2N 

D = e ' 0 

I = L = 1 


(27) 


a2 

where 1 min refers to the minimum distance among all unmerged differ- 
ence sequences of length N or larger, with h^^^ in effect at time n. 
Note the above is a symbol error probability and does not require 

differentiation as in the infinite memory case. This reflects the fact 
that only a single error is released when a truncation error event 

occurs . 

It remains to formulate the bound when the other K-1 deviations 
are in effect at time n. The matrices lQ (2 through are similar in 
form to that of (26), except the subscripts are cycled in straight-for- 
ward manner. Then similar bounds to that of (27) are found, each 

2 

using a different Adding all such truncation probabilities and 

dividing by K provides an upper bound on symbol error probability due 

p 

to truncation, ^^trunc' earlier argued the final upper-bound is 

obtained by combining the bound for infinte memory with Pg : 

° trunc 


P < P + P (23) 

® total ^ML ^ trunc 


where P is given by (25b). 
ML 
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V. NUMERICAL EXAMPLES ' 

The preceding development has been applied to the earlier exam- 
ple, as well as to two other codes: the binary {4/8, 5/8, 6/8} code 

and a 4-ary {3/16, 4/16} code. The free distance calculations for these 
designs project an asymptotic eiiergy efficiency relative to PSK or 
QPSK, of 1.4 dB, 2.8 dB, and -0.8 dB respectively. The first two 
have spectra comparable to that of CPFSK with h = 5/8, and the 99% 
bandwidth is approximately 1.67 times the bit rate. The 4-ary code hat’s 
a lesser bandwidth occupancy (the most rapid phase change is 3n/4 
radians per two bits, or that of binary signalling with h = 3/8. Its 99% 
bandwidth is approximately 0.7 times the bit rate, [8]. 

Figure 5 illustrates the bounds for the {2/4, 3/4} code, for unlimi- 
ted path memory and with N = 2 and N = 4, For unlimited memory, 

the upper and lower bounds closely bracket the true performance for 

-3 

typical levels of interest, say P < 10 . Simulation results for unlimi- 

ted memory actually lie close to the upper bound. All simulation results 
include at least 100 symbol errors. 

The decision depth (the depth in the trellis beyond which all 
unmerged paths have distance exceeding is four for this code. 

Thus we expect N = 4 should be quite adequate as a practical path 
memory length. Indeed the upper bound is only 0.2 dB displaced from 
the unlimited memory bound. Simulation results again confirm the 

validity of the bound. With N = 2, d^|^ is appreciably less, and in 
fact the efficiency is worse than that of PSK. This case indicates the 
effect of a poor choice of N. 
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Figure 6 presents upper-bound results for the K = 3 {4/0, 5/8, 
6/8} code. The larger asymptotic coding gain Is clearly revealed, and 
we also note this stronger code requires increased delay to realise its 
full potential. This is completely analogous to results for convolutional 
codes, where Increasing constraint length necessitates larger decoder 
memory. The decision depth in this case Is nine bits, and the upper 
bound for N ® 8 is asymptotically different than the unlimited memory 
bound. 

Finally, results for a narrow-band 4-ary design are shown In 
Figure 7. The decision depth in this case is = 9 s ymb ols . Again 
note that eiglit symbols of path memory is essentially adequate to 
achieve unlimited-memory performance. The asymptotic efficiency is 0.8 
dB poorer than PSK/QPSK; however the spectrum is quite compact. 
4-ary multi-li signalling appears to offer extra gain in tlie power/band- 
width/complexity trade-off, [9], [10], as is the case in other modula- 
tions such as CPFSK and partial-response FM. 
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3.0 4-ARY MULTI-H CODES 

By allowing non-binary signalling, we allow more freedom in the 
signal design and thus more flexibility in the power/bandwidth/complex- 
ity tradeoff. In more traditional designs, it is known that M-ary ortho- 
gonal signalling provides power savings over binary signalling, while 
M-ary PSK conserves in bandwidth. More recently, it has been found 
that non-binary partial response FM modulation affords substantial 
improvement in power-versus bandwidth performance, [1]. 4-ai'y and 
binary CPFSK are also known to be superior to binary CPFSK. 

We have surveyed 4-ary multi-h phase codes witli rectangular and 

raised-cosine frequency pulse shaping, [2]. It happens that to have 

the full constraint length of K + 1 symbols for a multi-h code with K 

deviations, that q, the common denominator of the rational deviation 

K 

indexes must be greater or equal to M , This is also the number of 
states in a maximum likelihood decoder. Consequently complexity grows 
rather rapidly with M, and we chose to stop at 4-ary, 2-h codes, where 
q = 16, 

Table 3.1 lists a selection of such codes having full constraint 
length, along with their free distance values. As for binary codes, the 
best codes are those with h's clustered near a common value, e.g. 
{3/16, 4/16}. Figure 3.1 illustrates the performance graphically as a 
function of h , the mean index. Note that for small fi, raised-cosine 
codes have slightly better distance (and slightly wider main lobe in the 
spectrum) and that a gain of about 2 dB over 4-ary CPFSK is typical. 
This corresponds to the typical gains shown for binary 2-h codes over 
binary CPFSK. 
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Table 3.1 Free distance of 4-ary 2-h codes with constraint length 3. 


TYPE 


MEAN 

LINEAR PHASE CODES 

R.C, FREQUENCY CODES 

Q 

'■i 

*■2 

INDEX 


72 

rree 

Decision 
Depth, Tb 


free 

Decision 
Depth, 1^3 

16 

4 

3 

.219 

/ 

1.688 

18 

/ 

2.103 

20 

16 

5 

4 

.281 

/ 

2.789 

24 


2.654 

22 

16 

7 

4 

.344 

/ 

3.454* 

30 

/ 

3.217* 

20 

17 

5 

3 

.235 


1.776 

22 


2.008 

20 

18 

8 

7 

.417 

/ 

4.527* 

34 


3.249* 

18 

18 

11 

8 

.528 


4.076 

12 


3.249 

10 

19 

8 

6 

.368 


4.262 

42 


2.560* 

28 

20 

4 

3 

.175 


1.105 

18 


1.387 

20 

20 

6 

5 

.275 


2.774 

34 


3.095 

36 

20 

8 

5 

.325 


3.399 

30 

/ 

3.549* 

42 

20 

9 

8 

.425 


4.798* 

44 


3.095 

12 

21 

10 

8 

.429 


4.286* 

52 


2.517* 

32 

22 

7 

6 

.295 


3.226 

46 


2.499 

34 

24 

5 

4 

.187 


1.323 

26 


1.633 

30 

24 

6 

5 

.229 


1.998 

36 


2.419 

42 

24 

7 

6 

.271 


2.774 

46 

/ 

3.217* 

52 

25 

6 

5 

.220 


1.853 

36 


2.249 

42 

25 

10 

7 

.340 


3.656* 

34 


3.726* 

40 

26 

5 

4 

.173 


1.136 

26 


1.406 

30 

27 

8 

7 

.278 


2.962 

60 


3.249 

62 

28 

8 

7 

.268 


2.780 

62 


3.285 

72 

28 

11 

9 

.357 


4.197 

88 


2.424 

52 

28 

12 

11 

.411 


4.646* 

76 


3.442* 

40 

29 

6 

5 

.190 


1.404 

36 


1.715 

42 

29 

13 

12 

.431 


4.838* 

80 


3.142* 

12 

30 

9 

8 

.283 


3.118 

78 ! 


3.095 

64 

30 

13 

12 

.417 


4.877* 

94 


3.586* 

22 

31 

6 

5 

.177 


1.238 

36 


1.515 

44 

32 

8 

7 

.234 


2.190 

62 


2.612 

74 

32 

10 

9 

.297 


3.423 

96 


2.654 

40 

32 

20 

19 

.609 


4.739* 

50 


3.579* 

36 


* : codes with best free distance anong codes with cannon q, 

/ ; codes of interest (having good power/bandwidth/coirplexitY trade-off) . 
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ipper bound 


A » / j 

2-h linear / 

2-h linear -J T ^ 


2-h R.C.\/7 , 

1 ! 



l-h linear 
'l-h R.C. 


o : 4 -ary 2-h linear phase cx>des 
A ; 4-ary 2-h raised-cxjsine codes 


Best free distance versus h. The codes are 4-ary 
codes with linear phase and raised-cosine pulse 
shaping. The unconnected marks represent codes with 
wide-spread indices (the difference between the numer 
ators of the indices is greater than one) . 
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The spectral properties of 4-ary 2-h codes may be calculated 
numerically using the procedure described earlier. An example is 
shown as Figure 3.2 for the {4/16, 5/16} case, Hsu [2] has found a 
quick rule of thumb for spectral behavior based on the mean index fi . 
Since the data values are (±1, ±3, ... ±M-1} and each symbol conveys 
log M 2 bits, we define an equivalent binary index as 

h. 


h. 


1 _ 1 _ 
K M/2 


K M/2 
1 1 (2i-l) 

j=l i=l 


log2M 


(3.1) 


= ( ) h 

'' log^M 

where R is the mean index, averaged over the K possible values. 
M-ary codes having an equivalent hj^ are found to have nearly equal 
spectral properties, i.e. 99% power bandwidth, etc. Since (M/2)/log2M 
= 1 for M = 4 and M = 2, we infer binary and 4-ary codes with equal fi 
will have comparable bandwidth. This is illustrated in Figure 3.3, 
showing Bgg and Bgg are equal within ± 10% when codes are compared 
on the basis of 

Table 3.2 summarizes bandwidth and distance statistics for a 
number of multi-h codes ranging from power efficient to bandwidth 
efficient designs. 

In summarizing the M-ary multi-h situation, we show performance 
in the distance-bandwidth plane. Figure 3.4. 4-ary 2-h designs show 
substantial gain over binary 2-h designs, particularly in the small 
bandwidth region, say Bgg Tj^ < 1. Also, we observe the rather de- 
pressing result that 4-ary single -h codes outperform binary 2-h codes 
under a bandwidth constraint. Since complexity is probably no worse 
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/ie,5/i6} 
/16,7/16} 
/16,9/16} 
/8,5/8,6/ 
16,5/16} 
31,10/31} 
/30, 13/30 
/32, 20/32 


Figure 3.3 90% and 99% bandwidtiis versus eqva valent binary mean 

index for various M-ary single index, linear phase codes. 
Also shown are some multi-h codes vdth ooirparable indices. 
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Table 3.2 Bandwidth performance for inulti-h codes with good 
power/bandwidth/cOTPplexity trade-offs. 


CODES 


"99 "90 "n» 


10.28 


100.5 


binary single h=5/7 
4-ary single lr=l/2 
binary {4/2,3} l.p. 
binary {4/2,3} r.c. 
binary {(5/3,4} l.p. 
binary {6/3,4} r.c. 
binary {8/4,5} l.p. 
binary {8/4,5} r.c. 
binary {9/4,6} l.p. 
binary {8/4, 5, 6} l.p 
binary {l0/5,6,7}r.c 

binary {8/2,3} r.c. 
binary {10/3,4} r.c. 
4-ary {16/3,4} l.p. 
4-ary {16/3,4} r.c. 

4-ary {16/4,5} l.p. 
4-ary {16/4,7} l.p. 
4-ary {16/4,7} r.c. 
4-ary {20/8,9} l.p. 


0.71 


0.78 


0.39 1.8 
0.43 1.8 


-0.38 


-0.74 


0.47 













for the 4-ary single-h design, it seems questionable whether binary 2-h 
codes, and even binary multi-h codes in general, should be studied 
further. We do note that 4-ary 2-h codes outperform 8-ary CPFSK in 
the small bandwidth region, and that the 4-ary 2-h class thus seems a 
strong modulation candidate. 
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4.0 CHANNEL SIMULATION 

We have implemented a FORTRAN software simulation of a typical 
satellite channel for purposes of assessing various effects on perfor- 
mance which are Induced by bandllmiting filters and/or saturating 
repeaters. Discrete-time simulation seems to be the expeditious ap- 
proach as analytical results of sufficient generality appear formidable. 
The simulation package is essentially completed at this point, though no 
error probability production runs have been completed. 

The system being addressed is shown in Figure 4.1. The software 
is highly-modular, so rearrangements of the functional blocks is possi- 
ble; however the configuration of Figure 4.1 represents the situation 
with post-modulation filtering and nonlinear amplification (either at the 
ground station or in a satellite transponder), down-link noise being 
dominant, and general receiver filtering. Coherent reception with 
maximum likelihood decoding is also assumed. In the absence of filter- 
ing and nonlinear distortion, performance is known for a variety of 
partial-response. FM and multi-h codes through calculations of minimum 
distance and error bounds, As with other more classical schemes 
however, the effect of chamml degradation is an important engineering 
question. In particular it is felt by some that these newer exotic 
designs are less tolerant of distortion than say QPSK or MSK. 

One may wonder why filtering is considered, since one of the 
attributes of continuous phase modulation is improved spectral shaping, 
with the. potential for avoiding the need for post-modulation filtering. 
It is however likely that, even so, some filtering will be performed for 
one of several reasons. First, equipment with formerly-necessary 
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channel filtering may still be in place and is to be used. Second, 
premodulation shaping alone is not able to attain spectral emission limits 
in high performance cases, e.g. trying to maximize data rate while 
meeting the FCC spectral mask. Post-modulation filtering to some 
degree is a simple expedient. Finally, aside from the above factors, 
the push toward high spectral efficiency will suggest moving bandwidth- 
efficient signals even closer in channel spacing, in essence placing even 
tighter specifications on spectral emission. 

The simulation processes the baseband complex envelope of a 
bandpass signal, and operates in discrete-time with eight samples per 
data bit , chosen to safely avoid aliasing phenomena. The signals of 
interest have spectra at least 30 dB down at the Nyquist frequency, 
which is four times the bit rate in this case. The data source is a 
255-bit maximal length shift register sequence, with an extra 0 ap- 
pended to balance the sequence. With eight samples per bit, we gener- 
ate 2048 complex samples as the signal sequence. These signal samples 
are used repeatedly with new noise sequences to generate enough trials 
to yield statistical significance in error probability estimates . This 
choice simulates much more rapidly than an approach which generates 
new signal samples continuously. 

The length of the PN sequence must be sufficiently long however 
to generate all significant unique data patterns as far as the channel 
output is concerned. Intersymbol interference is worse for certain 
patterns, and all significant patterns must be generated in the correct 
proportion. One of the properties of maximal length sequences is that 
the relative frequency of different N-tuples is very close to that of a 
coin-flipping process. With the degree of channel filtering expected for 
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our simulation, a L = 255 bit sequence is expected to be adequate, as 
the essential memory of the channel does not exceed 8 bits, and the 
relative frequencies are correct for up to 8-tuples. 

The receiver simulation adds white Gaussian noise to the received 
signal and filters the sum (if desired). The noisy discrete-time process 
is processed by complex correlators to yield statistics for the Viterbi 
algorithm receiver. The system as described has no carrier phase 
shift, so this need not be compensated. However, there is group delay 
roughly proportional to l/(composi.e filter bandwidth) which is on the 
order of one bit or more. This delay must be compensated by the 
correlators, and can be done once for the various configurations. Our 
timing quantization is Tj^/8. 

Once synchronized, the correlator outputs are passed to the ML 
routine having an appropriate state trellis for the modulation. After a 
sixteen bit delay, bit decisions are released by the decoder and com- 
pared with the actual PN sequence to compile error statistics. 

Note that the receiver is optimal for AWGN environments and acts 
as if no channel distortions were present. Equalization could undoubt- 
edly improve performance for severely bandlimited channels, and this 
could be investigated. However we feel that the initial study should be 
of the effect of channel distortion on performance of "conventional" 
receivers . 

4.1 SIMULATION RESULTS 

We first present results for MSK, a case for which results have 
recently appeared in the literature [1], [2], [3]. In Figure 4.2 we 
show the power spectrum at the modulator output, after filtering, and 
following hard limiting. The spectrum at the modulator output has side 
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lobes which decay as and well-defined nulls. The minor roughness 
in the spectrum is due to finite-averaging of a random process. The 
filter used was a 0.5 dB ripple, 4-pole Chebyshev characteristic with a 
BT product of 1.0. Note that this filtering truncates the mainlobe 
slightly, and obviously reduces sidelobe levels. The spectrum at this 
point would be comparable to QPSK with filtering at the first null in the 
spectrum (BT = 1). The interesting effect of ideal limiting is to re- 
generate spectral "sidelobes" having level almost as large as originally 
present. The power spectrum is smoother, and lobes have been era- 
dicated. This is completely in accord with recent results of [1] . 

Figure 4.3 presents similar results on spectra for a soft-limiter 
modelled as an error-function characteristic for AM-AM and a 3“/dB for 
AM-PM. This is perhaps more representative of typical TWT character- 
istics. Note with 3 dB backoff, the total power, relative to that of full 
output, is reduced, but more importantly the spectral regeneration is 
less pronounced. Figure 4.4 does the same comparison with a 9 dB 
backoff so that the characteristic is essentially linear; however the 
AM-to-PM conversion effect is still present. 

We also have capability for plotting amplitude and phase trajector- 
ies of the signal enroute. Figure 4.5 shows the amplitude, or enve- 
lope, of the filtered signal for filtered MSK with BT = 1. The bit 
sequence corresponds to 000000010111000. This plot gives an indication 
of the amplitude ripple produced by bandlimiting a constant-envelope 
signal, and in turn some notion of the severity of nonlinear amplifier 
distortion. Phase plots are of more interest as phase is the informa- 
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tion-bearing variable; however, the plots are currently difficult to 
interpret because of modulo-2n restrictions in the angle definition. 

Spectral results for double raised-cosine (DRC) with h * 1/4, 

channel filtering with BT = 0.5, and hard limiting are shown in Figure 

4.6. This design exhibits high bandwidth efficiency with Bgg T s 0.7, 

-8 

and sidelobes roll-off as f . We again note that limiting reconstitutes 
spectral sidelobes to a large degree. Figure 4.7 repeats the analysis 
for 9 dB back-off as above. 

A similar set of plots is shown in Figures 4.8 and 4.9 for 4-ary 
multi-h coding with deviation indices 4/16 and 5/16. The coding gain is 
1.4 dB and the ‘Bgg Tj^ = 0.9, 25% less than for MSK. 

The sidelobe reconstitution is an interesting phenomenon which has 
received considerable attention for QPSK and offset QPSK, Basically, 
when a signal is hard-limited, the rather smooth trajectories of the 
filter output are "sharpened", moving a relatively larger percentage of 
power to higher frequencies. (Only the phase trajectory varies at the 
output of the limiter). However the phase trajectories remain smoother 
than at the modulator output, so sidelobe levels are still reduced by 
filtering, then limiting. We have shown that the asymptotic rate of 
decay at the limiter output is equal that at the filter output; however, 
asymptotic rates tell little of absolute levels or where the asymptotic 
rate pertains. We expect to study this further under the next contract 
(NAG3-141). 

4.2 SOFTWARE DESCRIPTION 

The complete listing of the software mouw^.cj xs lounu in Appendix 
A. Program SIMULl simulates a digital communication system with 
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characteristics similar to those of a typical satellite link. It is com- 
prised of five basic sections; source (SORCEl, MARY), modulator 
(CPFSK, MULTIH, PARTLR), transmit filter (DIGBK), non-linear ampli- 
fier (TWT) and receive filter (DIGBK). Due to modularity, each sec- 
tion can be easily modified or totally skipped depending on the applica- 
tion. Furthermore, after each section the user can obtain complete 
documentation including envelope plot, and power spectrum phase tra- 
jectory, (PHASP, ENVELOP, FFTBK). When SIMULl is finished, 2048 
samples of phase, represented in complex form are written to a file for 
use withe the receiver program. The over-all flow diagram is shown in 
Figure 4.10. A description of each section follows. 

SORCEl and MARY 

Subroutines SORCEl and MARY comprise the data source for the 
simulation. SORCEl generates a length 256 psuedo-random sequence of 
I's and -I's and stores them in an array for future reference. The 
sequence is formed using a 8-bit shift register with feedback. Subrou- 
tine MARY converts the binary sequence generated by SORCEl into a 
M-ary sequence of length 256/K, where K = log 2 M and K is specified by 
the user. The M-ary symbol values are contained in the set 
{-(M-1), ... ,-3, -1,1,3, .. .(M-1)} and the mapping is arbitrary. 

CPFSK, MULTIH AND PARTLR 

Subroutines CPFSK, MULTIH and PARTLR form the modulator 
section of the simulation. In each case 2048 samples of the phase (8 
samples oer bit\ represented in complex form are returned in a 2048X2 
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GPFSK simulates a continuous-phase FSK modulator with deviation 
parameter h specified by the user. The psuedo-random M-ary sequence 
generated by the source, in conjunction with h, detei mines whether the 
phase should advance or retard and by how much. The phase varies 
linearly over 1 symbol time with accumulated phase m = dj^hll/M. For 
the special case h = .5 the modulation becomes MSK. 

MULTIH simulates a multi-h phase code modulator with hj^,h2...hj^ 
specified by the user. MULTIH is similar to GPFSK except that h is 
now time varying. 

PARLTR simulates a partial response FM code modulator with the 
response time K, the deviation parameter h and tl\e pulse shape g(t) 
specified by the user. The pulse shape g(x) can be either rectangular 
or raised-cosine and lasts a total of K symbol durations. The composite 
phase at any one sample point is thus a function of the phase induced 
by the present symbol d^^ and the phase Induced by K-1 symbols of the 
past. For the special case K » 2 and gCx' rectangular (binary source), 
the modulation is known as duobinary FM. 

DI GBK 

Subroutine DIGliK simulates a low pass, 4-pole Ghobychev digital 
filter with bandwidth BW specified by the user. The filter coefficients 
arc prc-calculated by a digital filter design program employing the 
bilinear transformation and then stored on disk for future use. Sub- 
routine DIGFIL within DIGBK reads these coefficients off the disk and 
filters the time series data via a recursive difference equation tech- 
nique. The time series data is run through the filter twice to establish 
steady state conditions. 
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TWT 

Subroutine TWT simulates a non-linear amplifier with characteris- 
tics similar to those of a Travelling Wave Tube amplifier used in may 
satellite transponders. Output distortion is a combination of two ef- 
fects, AM- AM conversion and AM- PM conversion. 

With AM- AM distortion, the output envelope is a non-linear func- 
tion of the input envelope. The non-linearity is either a hard limiter in 
which case the output envelope is set to unity regardless of the shape 
of the input envelope, or a soft limiter described by the equation y = 
erf(x/a). In the latter case a controls the degree of limiting; i.e. as a 
increases the amplifier becomes more linear and as a decreases the 
amplifier saturates and behaves more like a hard limiter. Erf(x/a) is 
computed using a 5-term power series expansion with maximum error 
less than 7.5 x 10“^. 

In AM-PM conversion, fluctuations in the input envelope cause 
phase distortion in the output waveform. This distortion is modeled by 
a 3° per dB characteristic, i.e. 

0 = -(3.0/57.3) 20 log(A) 

where O is the resulting change in output phase and A is the magnitude 
of the input waveform. 

FFTBK, PHASP and ENVELP 

Subroutines FFTBK, PHASP and ENVELOP comprise the documenta- 
tion section of the simulation. FFTBK computes and plots an estimate 
of the power spectral density, S(f), while PHASP computes and plots a 
phase trajectory (modulo 2n) and ENVELOP computes and plots the 
envelope or magnitude of the complex time series data. 


IS 
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Within subroutine FFTBK the time series data is first pre-multi- 
plied by a Hanning window of the form 

H(n) = .54 + .46 cos C|^|g) 

to improve accuracy at low values of S(f) (-50dB down). The modified 
time series data is transformed into the frequency domain via a 2048 
point complex FFT. The magnitude of the resulting, frequency samples 
are then averaged with 10 adjacent samples to provide a smoother 
estimation of the power spectral density. 
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C 

C 

c 

c 

c 

c 


c 

c 

c 

20 


212 


i:l4 


C 


198 


PROGRAM SIMIJLI 


THIS PROGRAM SIMULATES A DIGITAL COMM» SYSTEM 
INCLUDING THE FOLLOWING COMPONENTSP SOURCE(DINARY OR M-ARY) 
MODULATOR ( CPFSK !> MULTI-H y PARTLR ) c XMTR FILTER POLE CHEDYCHEV 
TWT NONLINEARITY (WARD OR SOFT LIMITER). RCVR FILTER- 


REAL MOD ID (4) rW (3) 

DIMENSION A ( 2048 . 2 ) . IBIT ( 256 ) . AR ( 2048 ) . AI ( 2048 ) 
EQUIVALENCE ( A ( 1 y 1 ) . AR ( 1 ) ) y < A ( 1 y 2 ) y AI ( 1 ) ) 

BYTE ANS y YEA y TMPFIL ( IS ) y CORFIL ( .1. S ) 

BYTE TFNAME ( 15 ) y RENAME < 15 ) y FNAME ( :l, 5 ) y TWTID ( 4 ) 
DATA FS/S ♦ 0/ y NDATA/2S6/ y IEXP/1 1 / y MRY/1/ 

DATA FNAME /'N' y '0" y "N' y 'E^ y 11>K' V 
DATA TFNAME /'N' y '0' y 'N' y 'E' y ll;i<' V 
DATA RENAME / 'N ' y ' 0 ' y ' N ' y ' E ' y 11-Y ' V 
DATA TWTID /'N' y '0' y 'N' y 'EV 
DATA YEA/'YVyN/0/yW/3>!<0.0/yNT /()/ 

DATA TMPFIL /'T' y "M" y 'P' y 'F' y M' y "L" .9^' V 
DATA CORFIL /'C ' y '0 ' y "R ' y 'F ' y ' I ' y ' L ^ y 9>r V 
DATA MODID /At' '/ 

CALL ASSIGN (50y TMPFIL y 14 y ' NEW ' ) 
ND-IFIX(FS)miATA 


GENERATE 255 BIT PN BINARY OR M-ARY SOURCE 


CALL SORCEl ( NDATA y IBIT ) 

WRITE( 7.212) 

FORMAT ( 'OBINARY OR M-ARY? (0 OR 1)') 
READ(5y214) IANS 
FORMAT (15) 

IF ( IANS ♦ EQ » 1 ) CALL MARY ( NDATA y IBIT y MRY ) 
CHOOSE MODULATION TYPE 
WRITE (7.198) 

FORMAT ( ' OCHOOSE MODULATION ' / 

TYPE'/ 

&' If CP-FSKV 

Z' 2f MULTI-HV 

8' 3f PARTIAL RESPONSE'/ 

S'$WHICH?' ) 
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READ (5y200> XSW 
FORMAT (16) 

lT(ISW,Eat.1.) CALL CPFSK(AyNi:iyNDATAyFSyIBITyl-l(X) yN 
8rMODKiyMRY) 

XFdSW .ECU 2) CALL MULTXH(AyNi;iyNi:iATA yFSy XBXTyHyNyMODXD 
8y|<yMRY) 

XF(XSW*E0,3) CALL PARTLR<AyNDyNDATAyFSy XBXTyl-Kl ) yN 
SyMODXDyNTyMRY) 

OTHER BRANCHES MADE HERE AS THEY BECOME AVAXLABLE 
IF(N ,NE* 0) 00 TO 27 

CALL PRXNT('NOT IMPLEMENTED AT THXS TIME') 

GO TO 25 

DOCUMENTATXON (PHASE TRAJECTORyFFTyENVELOPE PLOT) 

WRXTE(7y220> 

FORMAT ( 'OPHASE TRAJECTOR PLOTTO 
READ (5y 204) ANS 

XF(ANS*EQ*YEA) CALL PHASP( A y AR »• AX yND) 

CALL FFTBK ( A y AR y AX y ND y XEXP ) 

XMTR FILTER AND DOCUMENTATION 

CALL PRINT ('DO YOU WANT CHANNEL FILTERING AT THE XMTR? (Y/N)') 
READ (5y204) ANS 
FORMAT (A:l.) 

IF (ANS .NE. YEA) 60 TO 30 
CALL D I GBK ( A y AR y A I y ND y TFNAME ) 

WRITE (7y 206) 

FORMAT ('ODO YOU WANT AN ENVELOPE PLOT (Y/N)') 

READ (5y204) ANS 

IF (ANS .E0» YEA) CALL ENUELP ( A y ARy AX y ND y FS ) 

WRITE(7y220) 

READ(5y204) ANS 

I F ( ANS ♦ EQ , YEA ) CALL. PHASP ( A y AR y A I y ND ) 

CALL FFTBK ( A y AR y AX y ND y XEXP ) 

TUT NOLXN AND DOCUMENTATION 

C A L. L P R' I N T ( •' D 0 Y 0 U W A N T A T W T' S I M U L A T' I ON? ( Y / N ) ' ) 

READ (5y204) ANS 

IF (ANS ♦NE. YEA) GO TO 50 

CALL TWT(AyNDyTWTID) 

WRITE (7y 206) 

READ (5y 204) ANS 

IF (ANS .EG. VEA) CALL ENUELP ( A y AR y AX y ND ) 

CALL FFTBK ( A y AR y AX y ND y XEXP ) 

RCOR FILTER AND DOCUMENTAION 
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SO CALL PRINK 'WOULD YOU LIKE ADDITIONAL CHANNEL FILTERING <Y/N)') 

READ (Sr 204) ANS 
IF (ANS «NE» YEA) GO TO 60 
CALL D I GBK ( A r AR r A I r ND r RFNAME ) 

WRITE(7r206) 

READ (Sr 204) ANS 

IF (ANS ♦EO. YEA) CALL ENUELP ( A r ARr AI r ND) 

WRITE (7r 220) 

READ (Sr 204) ANS 

IF(ANS»EO»YEA) CALL PWASP(Ar ARr AI rND) 

CALL FFTBK ( A r AR r AI r ND r I EXP ) 

C 

C WRITE TIME SERIES DATA TO A FILE 

C 

CALL CLOSE (50) 

60 CALL PRINT ('DO YOU WANT TO WRITE SAMPLES TO A FILE? (Y/N)') 

READ (Sr 204) ANS 

IF (ANS .NE, YEA) GO TO 70 

CALL PRINT ('INPUT FILENAME') 

IERR-0 

CALL GETSf R ( 5 r FNAME r 14 r lERR ) 

CALL ASS I GN ( :l 0 r FNAME r 1 4 r ' NEW ' ) 

WRITE (10 ) FNAME r MODID r NT r H r FS r TFNAME r TWTID r RFNAME r K r IBIT r A 
CALL CLOSE (10) 

70 WRITE(7r203) FNAMEr (MODID( J) r J=-"l r2) rNTr (MODID( J) r J=3r4) 

WRITE ( 7 r 209 ) H ( 1 ) r H ( 2 ) r H ( 3 ) r FS r TFNAME r TWTID y RFNAME 

208 FORMAT ( 'OSUMMARY OF SIMULATION'/ 

X'OFILENAME: 'rlSAl/ 

%' MODULATION TYPE: ' r2A4r3Xr I2r2A4) 

209 FORMAT( ' W( 1 )==••'»• F8 ♦ 3r 3Xr 'W(2)--' yF8»3r3Xr 'l-l(3) = ' rF8*3/ 

X' FS~'rF8»3/ 

8' XMTR FLTR: 'rlSAl/ 

X' TWT : 'v4Al/ 

8' RCVR FLTR: ' rlSAl) 

WRITE(7r210) 

210 FORMAT ( ' OSIMULATION IS COMPLETE') 

STOP 

END 
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C 

C 

C 

C 


C 

C 


30 


50 




SUBROUTINE SORCEl (NBATAr IBIT) 

THIS PROGRAM SIMULATES A 255 BIT PN GENERATOR 
NBATA=>i OF SOURCE BITS? BIT- BATA? Cl«d.?0“ -I) 


DIMENSION IBIT<Ni:iATA) 
LOGICAL D<8)?F(3> 

DATA IU/32757/ ? XX/7057/ 
DATA D/7>!< ♦ FALSE ♦ ? . TRUE » / 
DO 50 I"=lrNDATA 
F(l)™( .NOT.D(S) ♦AND.D<7) ) 
F(2).”( .NOT.Fa) t AND,ri(3) ) 
F(3>'«( .N0T.F(2> .AND,D(2) ) 


♦OR* (D<S).AND.(«N0T.D(7)>) 
♦ OR* (F(l) ♦AKD« < ♦N0T»D(3) ) ) 
♦OR* <F<2) *AND* ( *N0T*D(2) > ) 


PERFORM THE RIGHT SHIFT 
DO 30 J™:L?7 
D(9"J)-D(8“J) 

CONTINUE 

D(l)?-F(3) 

IF (D(8) *E0* *FALSE* ) IBIT(I) 

IF (DCS) *E0* * TRUE* ) IBIT(I) - I 
CONTINUE . 

RETURN 

END 
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SUBROUTINE MARY < NDATA y IBIT r M ) 

C 

C THIS SUBROUTINE MAPS BINARY DATA INTO M--ARY DATA 

0 

DIMENSION IBIT ( NIJATA ) 

20 WRITE(7y200) 

200 FORMAT ( 'OINPUT K ( M""2:K?|!K ?M AS IN M-ARY>V 

X' 1< K >- 6 ') 

READ <Sf2I0) M 
210 FORMAT (IG) 

IF<MJ.E,1) GO TO 20 

IF(M,GE.6^ GO TO 20 

C 

C PERFORM MAPPING 

C 

L~0 

K“0 

MINDIT=~(2>Ka'<M)KL 
30 DO 80 K -If NDATA yM 

L»LT1 
J^MINDIT 
DO 60 li^'^lvM 

IF ( IDIT ( lvM™l ) ♦ EQ . ) IDIT ( KTI-l )= =0 

JT ( 2m ( M-”ITI ) ) >KIDIT < KTI-;I. ) 

60 CONTINUE 

IBIT(L)™J 
SO CONTINUE 

RETURN 
END 
5K 
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C 

c 

c 

c 

c 

c 


c 

c 

c 

100 


102 


c 

c 

c 


10 

20 


SO 


SUBROUTINE CPFiSK ( A . NB . NBATA r FS y IBIT y H y N r MODID y M > 

=.rrs: sss 's,r =s.ss-:, 

SSft SS U SSS S S f wSSSS ” ““ 


DIMENSION A(Ni:iy2) y IBIT(NDATA) yTODTIU?) 
REAL PIyM0DID(2) .'.i. w.; 

DATA TODID/ ' CP-F ' y ' SK ' / 

DAT A PI/3 , 1 4 :l. 592/ y ANf3/0 ♦ / y K/ 0/ 

DATA UFST/0»0/ 

MODIDa).'=TODIDa) 

M 0 DID( 2 )™'r 0 DID( 2 ) 


DETERM INE MODULATOR PARAMETERS 
WRITE (7y 100) 

FORMAT ('OTHIS SIMULATES CP-FSK. V 
$' INPUT H') 

READ (5yl02) H 
FORMAT <F6.3) 

IFS.-=IFIX(FS);f'M 

NSYMBL--NDATA/M 


P E R F 0 R' M S I M ( J L A 7‘ 1 0 N 

I-O 

l-X-i-l 

J ™0 

j~j-n 

ANO “OFST'f ( I BIT ( I ) >KH:KPXh^J > / ( FS>!<M ) 
K“K +1 

A ( K y 1 ) :.-00S ( ANO ) 

A ( Ky 2 )«SIN(ANO) 

IF(J.LTMFS) 00 TO 20 
OFST^.'nOFSTMBIT( I ) :{<H:i:PI 
IF(IJ..T,NSYHBL) 00 TO 10 
IF(K»E0,ND) 00 TO 50 
IFS=!=ND-”-IFS:{<NSYMi;il.. 

00 TO 10 
CONTINUE 
N=".l 

RETURN 


6'4 


7 


p/P/f rvt, , 


QF f'Of ” 




c 

c 

c 


c 


c 

c 

c 

w 

200 


':>7 

202 


204 

30 

206 

208 

C 

C 


50 


90 

."ic 


SUBKOiiTINE MULTIH ( A y NB y NBATA y FS y XBITy H y N y MOBIB y Ky M ) 

ry:™ 

BlMENblON A<NBy2) y|-l(3> ? I’BIT (NBATA) yTOBTB(''>> 

REAL M0BIB(2) " y y , uuxKU, ;> 

MOBIBCDn-'TOBIBCl) 

MOBIB (2).~tob:i;b(2> 

N«.'=l 

NSYMBL-NBATA/M 

IFS«IFIX(FS)>f;M 

I NPUT PARAMETERS ( K y H ( K ) ) 

WRITE(7y200) 

I ' ZZr^t ho«ul„torv 

8' (NOTEy KMAX 3)') ~ ^ 

REAB(5y202) K 
rORMAT(I3) 

1F(K,LT.4) go to 30 
WRITE(7y204) 

™RMAT( '0KMAX «: 3y INPUT K' ) 

WRITE(7y206) 

FORMAT ( ' OINPUT H (I ) y , * . H f f; ) ' ) 

R L A B ( 5 y 2 0 8 ) ( |-| ( I ) y l .•=.• ;|. ,, > 
format (3F6. 3) 

PERFORM SIMULATION 

OFST-'^O.O 
1^=0 
I™I + ;l. 

MH:<=M0Baf-(K~..;|,),|Of;l. 

J™0 



ANG«=OKSTf ( IBIT ( I ) ;M~i ( MH ) >fcpl:,<< j ) / ( ^ 

L“"L 'M. 

A ( L y I ) :«C0S ( ANG ) 

Aa.y2)«»-SIN(ANG) 

IF(J J..T.IFS) GO TO 75 

0 1" S T "•••• 0 F S T •}• I B I r < I ) ;K |-| ( m H ) $ p i 

I F ( I , L T ♦ N S Y M B L ) G 0 T 0 0 

IFCL.EO.NB) GO TO 90 

I FS=NB- 1 FSi'cNS YMBL 

GO TO 50 

RETURN 

ENB 
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C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

n 


200 


30 

0()0 
Am 1/ A4> 


•'-iO 

20'1 


SO 

208 

2 : 1.0 


SUDROUTINE PhRTLR < A r Ni;i y NDATA y FB y XB-U y H y N y MOnXB y NTT y MRY > 


THX8 RR08RAM 8XMULATE8 A PARTXAI. REBPONBli: FM COBE 
MODULATOR FOR TWO TYPE8 OF PUL.BEOy DUOBXNARY FM 
(N-RECTANOULAR) AND DOUBLE RAXOED C08XNE (NRC) 

A™TXME SERXE8 BATAy ND- # OF DATA P0XNT8.- NDATA-# OF 
DATA BXTBy FS-SAMPLXND RATEy XBXT -DATA BXT ARRAY y H~ 
DEOXATXON PARAMETER y N XS BET TO iy MODXD^ XDEN* 
PONSE TXMEy MRY AB XN M-ARY 


NTT--REB" 


D X MEN S X 0 N A < N D y 2 ) y X B X T ( N D A T A > y T 0 D .X D ( <f> ) y X R I" M < 7 ) 


REAL M0DXD(4) 

COMMON PI'S y F* X y PM v NT y MM 

DATA TODXD / " PART ' y ' REBP ' y ' -REC) ' 

DATA P X /3 * 1 4 :l. SY26S3S 9/ y K/O/ 

DATA NTMAX /?/ 


•RC ' y ' 


PFS»“FB 

N-;l 

MODXD<X)=sTODXDa) 

M0DXD(2>-T0DXD(2) 

MM-MRY 


V 


DETERMXNE MODULATOR PARAMETERS (PULSE 'SHAPE ylO 


WRXTE (7y 200) 

FORMAT ( 'OTMXS SXMULATES A PARTXAI... RESPONSE FM MODULATOR'/ 
X' FOR TI!E FOLLOWXNS CASES'/ 

X' :l. . DUOBXNARY FM < N--RECTANOULAR > ' / 

X' 2v DOUBLE RAXSED COSXNE (NRC)'/ 

X' WHXCM?') 

READ (Sy 202) XSW 
FORMAT (X6) 

XF(XSW.GT«2) 00 TO 40 
XF(XSW*Lr*:l.) so to 40 


(30 TO SO 
WRXTE(7y2()4) 

FORMAT ( 'ONO SU(.:H CI-IOXCE y 


CHOOSE AOAXN' ) 


00 TO 30 
WRXTE ( 7 y 208 ) 

FORMAT ( 'OXNPUT HyNT'/ 

X' (NT-T'iURAllON OF FRECn PULSEy NTMAX™/)') 
READ (Sy 210) HyNT 
FORMAT ( FA *3yXS) 


MTT™NT 
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IF<NT.LT*NTMAX) (3(3 TO (30 
WRITE (7 I'll. 9 (i)) NTMAX 
198 FORMAT <" NT MIJ(3T BE <SIS) 

(30 TO 00 
80 CONTINUE 

IF<ISWJ-Q*2) f30 TO 80 

modid(3)-tod:i:i.uo) 

M0tiILU4)”T0i:iKi(4) 

GO TO 87 

80 M0rinK3).~M0D;i:ii(0) 

M0BID(4)«-'T0i;ini(6) 

87 CONTINUE 

PH--W 

IFS==^IFIX<FS)-'f(MM 

NSYMBL.™NDATA/MM 

C 

C PERFORM SIMULATION 

C 

BO 90 IN«":UNT-;l. 
IPREV<NTT:l.-IN)™:l. 

90 CONTINUE 

IPREOa)-<i;BITCI. ) 

0^=0 ♦ 0 
M="0 

3.00 M=-M-fl 

J~0 

no j™jf3. 

K-KTl 

L-J 

P=.-0*0 

BO 120 IN“3, i»NT 

P"P-}-PHASE ( IPREV ( IN ) y L j- ISW ) 

L“LT(8;}:MM) 

3.20 CONTINUE 

P--BTP 

A < K n ) ™COS ( P ) 

A(Ki-2)“(3IN(P) 

IF(JJ_T. IF(3) (30 TO 3.3.0 
B •••- B •{• p i-i ;k F' ;i; 1 1"' e o ( n T) 

BO 135 IN-lyNT-1 
IPREO ( NTTl -IN ) «=IPREO ( NT-IN ) 
130 CONn.>UE 

IF(M«C)EtN(3YMBL) (30 TO 140 
IPREO(n-=IBIT(Mfl) 

GO TO 100 

140 IF(K.E(,3»NB) GO TO 100 

IFS™NB™IF(3.^NBYMBL. 

IPREVa )™IBITCM-M ) 

GO TO 100 
150 RETURN 

RETURN 
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END " 

REAL FUNCT I ON PHASE ( L t- M y N > 

COMMON PFSi^PIrPI-UNTyMM 
iF(Nj=:a*:i.) (so to so 
PHA» < PHJKPIJKL ) / ( NTiKMM ) . 

phas-“S:l'n ( ( 2 . o!i<p:i:>i<M ) / ( pfs«nt>kmm > > 

PHASE«PHAH< ( < M/PFS )•-■((< NT>KMM ) / ( 2 ♦ 0:{<P]: > ) JKPHAS ) ) 
RETURN 

PHASE« ( L.?!<PH>!<P:DltM ) / ( NT:KPFS5«MM ) 

RETURN 

END 
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SUBROUT I NH- ENUELP ( A s- AR y A 1 * ND y FS ) 


THIS PROORAH PLOTS A GRAPH OF THE SIGNAL 
ENVELOPE VS* TIME, A>"TIME SERIES DATAy ND= 

# OF DATA POINTS y FS==SAMPLING RATEy AR(K>»A(Ky J. ) 
AI(K>~A(K'y2) 


DIMENSION A < ND y 2 ) y AR ( ND ) y AI ( ND ) 
REAL-KS XNAME(2) yYNAME(2) 

DATA XNAME /' TIME 'y' SECONDS'/ 

DATA YNAME/ ' S I GN AL ' y ' ENVELOPE ' / 

NS=100 

NSMAX=150 

WRITE (50) A 


DETERMINE PLOT PARAMETERS 
WRITE (7y 200) 

FORMAT ( 'OACCEPT (<CR>) OR CHANGE (/<NEW VALUEXCR 
8' NUMBER OF SAMPLES TO BE PLOTTEIKNS) V 
X' ( REMEMBER y # OF PAGES™NS/5.1. ) ' ) 

CALL I NT 1 0 < ' NS ' y NS y NSMAX y 0 ) 

CALCULATE MAGNITUDE 

DO 10 J-:l, yNS 

AI ( J ) -SORT < A < J y 1 > «<^<2fA ( J y 2 ) >lo!<2 ) 

AR ( J ) =•• J/FS 
CONTINUE 

PLOT ENVELOPE ON CRT OR DECWRITER AND SAVE DATA 
ON FILE FOR HP-PLOTTER 

PAUSE 'ROLL PAPER TO NEXT TO BOTTOM LINE THEN <CR 
DUMB™0,0 

CALL PLOTR ( AR y A I y NS y XNAME y YNAME y 0 . 0 y 1 . 25 y 1 ) 

CALL HPPLT(AyNDyNS) 

REWIND 50 
READ (50) A 
REWIND 50 
RETURN 
END 
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C 

c 

c 

c 

c 

c 


c 

c 

c 


200 


202 

C 

C 

C 


50 

C 

C 

C 

C 


SUBROUTINE PHASR ( A i- AR y AI i- ND ) 


THIS SUBROUTINE CALCULATES ANB PLOTS PHASE TRAJECTORIE 
PHASE :== ARCTAN(AI/AR) MODULO (2PX) 

A“TIME SERIES DATAy ND«# OF DATA POINTSc AR(K)==A(Ky 1 ) 
AI(K)“A(Ky2) ♦ 


DIMENSION A ( ND >■ 2 ) v AR ( ND ) y AI < ND ) 

REAUK8 XNAME < 2 ) y YNAME ( 2 ) 

DATA PI/3 ♦ 14:1.592654/ 

DATA XNAME/ ' PHASE ' v ' ( RAD ) ' / 

DATA YNAME/ 'SAMPLE 'y 'INDEX '/ 

WRITE (50) A 

DETERMINE PLOTTING PARAMETERS 

TP=2!KPI 
IFP™1 
NP-150 
WRITE (7y 200) 

F 0 R M A T ( ' 0 F 1-1 A S E T A J E C T' 0 R Y P L 0 T T E R ' / 

&' ACCEPT (<CR.>) OR CHANGE (</NEW UALUE><CR.>) '/ 

S' THE # OF POINTS TO BE PLOTTED (NP') 

CALL INTIO ( ' NP ' y NP y ND y 0 ) 

WRITE(7y202) 

FORMAT( 'OFIRST POINT TO BE PLOTTED ( IFP) ' ) 

CALL INTIOC ' IFP' y IFPyND-NPyO) 

CAALCULATE PHASE MOD (2 PI) 

DO 50 J::>:IFPyIFPTNP 
AI(J)==ATAN2(AI(J) yAR( J) ) 

AR(J)™J-1 

CONTINUE 

PLOT PHASE TRAJECTORY ON CRT OR DECWRITER AND SAVE DATA 
ON FILE FOR USE WITH HP-PLOTTER ( SUBROUTIN HPPLT) 

CALL HPPLT(AyNDyNP) 

PAUSE 'ROLL PAPER' 

CALL PLOTR ( AR y AX y NP y XNAME y YNAME y "-4>KPI y 4-'!<PI y 0 ) 

REWIND 50 
READ (50) A 
REWIND 50 
RETURN 
END 
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C 

c 

c 

c 

c 

c 


c 


202 

204 


C 

C 

C 


20 

C 

C 

C 


206 


SUBROUTINE FFTBK < A y AR i- A 1 9 ND y I EXP > 

THIS SUBROUTINE COMPUTES AND PLOTS A POWER SPECTRAL 
BENSITY VIA A HANNINO WINDOW FOLLOWED BY A FFT AND A 
N--POINT SMOOTHER* A«TIME SERIES DATAy ND“# OF BATA POINTS 
IEXP -4-0G2 ( ND ) y AR ( K ) ™A ( K y I ) y AI ( K ) ( K y 2 ) 

REALMS XNAME(2) yYNAME(2) 

DATA XNAME/' (F/R) "y' V 

DATA YNAME/' GAIN 'y' (DB) V 
DATA PI/3* ;14;IS9265;JSS/ 

DIMENSION XCJ.29) 

DIMENSION A ( ND y 2 ) y AR ( ND ) y AI ( ND > y B ( 2 ) y C ( 2 ) 

BYTE ANSyYEA 
DATA YEA/'YV 


WRITE(7y202) 

FORMAT <'0D0 YOU WANT A FFT? (Y/N>'> 

READ <Sy204> ANS 
FORMAT (A I) 

IF (ANS *NE* YEA) RETURN 
WRITE (SO) A 

MULTIPLY TIME SERIES DATA BY HANNIN(.v WINDOW 
DO 20 K-'IyND 

Y=FLOAT ( K-ND/2 ) /FLOAT ( ND ) 

W«0 » 5T0 * 5.YC0S ( 2:{JP I :K Y ) 

A(Ky.l)-~A(Kyl )!KW 
A(Ky2)==A(Ky2)^W 
CONTINUE 

CALL FFTR ( A y ND y I EXP y ~ 1 * ) 

DETERMINE SMOOTHING AND PLOTTING PARAMETERS 


KS=11 

NLP=5I2 

NND=---ND/2 

NBP==4 


WRITE(7y206) 

FORMAT ('OACCEPT ( <CR> ) OR CHANGE (/<NEW OALUEXCR> THEV 
FOLLOWING PARAMETER UALUES ASSOCIATED WITH THE PLOT'/ 
X' # OF SAMPLES BETWEEN PLOT POINTS (NBP)J') 
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208 


C 

C 

C 


os 

XL vJ 


26 


28 

C 

C 

C 

C 


t 


CALL IMTIOC' NBP' »NBP»40>0) 

CALL PlJTSTR(7i'' LAST SAMPLE POXNT TO BE PLOTTED ( NLP ) J 0 ) 
CALL INTIOt' NLP' yNLPyNNDvO) 

WRITE ( 7 y 208) 

FORMAT ( ' OSMOOTHINO CONSTANT ( KS ) ' / 

X' KS MUST BE ODDyKS”! FOR NO SMOOTHING') 

CALL INTIO<' KS' rKSrlSyO) 

KSl = (KS“-l)/2 
KS2=:KS1 + 1 
NPTS=NLP/NBP 


CALCULATE MAGNITUDE AND PERFORM KS-fOINT SMOOTHING OPERATION 


DO 25 J=:l.rNND 

AR ( J ) =A ( J »• ;L ) :K>K2+A ( J i- 2 ) ;{c$2 

K=0 

DO 27 J==:-l»NLPyNBP 
K“K+1 

X(K)“04 

IM=--“MAXO( J-KS;U'I) 

DO 26 KK~IM^' J-f-KSI 
X(K)~X(K)fAR(KK) 

IM2-MIN0 ( 5 1- J-M<S2“-IM ) 

X < K ) "-X ( K ) /FLOAT ( IM2 ) 
x<k)“:lo.^alog:i.o(X(K) ) 

IF ( X ( K ) ♦ LT , -50 ♦ ) X ( K ) “-50 » 0 
AR ( K ) ==FLOAT ( J-1 ) /256 ♦ 

DO 28 KK=:l. ^NPTS 
AI(KK)-"=:X(KK) 


PLOT POWER SPECTRUM ON CRT OR DECWRITER 
FILE FOR HP-PLOTTER* 


AND 


SAUE DATA ON 


CALL HPPLT(AyNDyNPTS) 

PAUSE 'ROLL PAPER TO NEXT TO BOTTOM LINE THEN <CR>' 
DUMB=0,0 

CALL PLOTR ( AR y AI y NPTS y XNAME y YNAME y -SO . 0 y 50 ♦ 0 y I ) 

REWIND SO 

READ ( 50 ) A 

REWIND 50 

RETURN 

END 


o o n n o o n 
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5 

C 

C 

C 


50 

200 


210 



230 

C 

C 

C 

C 

70 


90 


SUBROUTINE TWT ( X y ND » TWTID ) 


THIS SUBROUTINE SIMULATES A pntMTr 

TYPE AMPLIFIER* X™ TIME SERIES DATAy NB=--#Or- BATA POINTS 

TWTIB^^IBENTIFICATION 

am-am determined by function 1-1 
AM-PM determined by function F2 


DIMENSION X(NDy2) 

BYTE TWTID(4) yDWTID(4) 

DATA DUTID /' Y' y 'E' y ■ 3' y ' V 
DATA S/1*0/ 

data PI/3* 14159265359/ 

COMMON SySP 
SP=1*0/SQRT(2>!!PI) 

DO 5 J"ly4 
TWTID( J)“DWTID( J) 

CONTINUE 


DETERMINE LIMITER TYPE <AND> SK3MA 


WRITE(7y200) 
FORMAT ( 'OHARD 


,I MITER OR SOFT I.IMITER? 


READ (5y 210) ISN 
FORMAT (41) 

IF(ISW*E0*1) GO TO 60 
IF(ISW.NE*0) GO TO 50 


GO TO 70 
WRITE(7y220) 

FORMAT ( 'OINPUT SIGMA' ) 
READ (5/230) S 
FORMAT ( F6 * 3 ) 


(0 OR 1) 


) 


P E R F 0 R M S I M U I... A T 1 0 N 

ISW---1 IMPLIES SOFT LIMITER 

DO 100 I=-=lyND 

ASSORT ( X ( I y 1 ) ;i<;{<2TX ( I y 2 ) ) 

PH I ■•y~ ATAN2 ( X ( I y 2 ) y X ( I y 1 ) ) 
AMP^A/S 

A0“F1 ( AMI"' y ISW ) 
DPHI==F2(AyISW) 

X ( I y 1 ) -AO>KCOS ( PH IP DP HI ) 
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X ( I » 2 ) -•■••AO>KSl’N < PHIMiPHJ; ) 

CONTINUE 

RETURN 

END 

REAL FUNCTION Fl<A>-IOW) 

IDEAL LIMITER 

DATA B.I. /♦3;l.9;:!:ei530/s- D2/“, 356563702/ 

DATA B3/:l. . 701477937/ x B4/™;l. ♦ 02:1.255970/ x B5/1 * 330274429/ 
DATA P/ ♦ 231641.9/ x PI/3 ♦ 1.4.1.59265359/ 

COMMON Ox OP 
IFdSW.ECnO) GO TO 50 
z~sp>KEXP ( " ( mrz> /2 . o ) 

Tl-I,0/(l.*0fP;i<A) 

T2-"T1>K-K2 

T3“T2«T1. 

T4-T35KT1. 

T5-T4>KTI 

F==2>K ( BI >}:T.I. •f-T2:KB2TB3-'KT3+B4>{{T4TB5:{<T5 ) 

F1.“2>K( »5~F) 

RETURN 

F1.™:U0 

RETURN 

END 

FUNCTION F2(AxISW) 

IF<ISW.EQ,;l. ) GO TO 50 

F2“- ( 3 ♦ 0/57 » 3 ) t ( 20;RALOG ;I. 0 ( A ) ) 

RETURN 

F2™0»0 

RETURN 



c 

c 


c 

c 


c 

c 

c 


200 


208 


:lo 


20 

C 

C 


C 


C 


30 
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SUBR8UTI NIi- D l’(3BK < A x AR y A I y ND y r-NAME ) 


THIS SUBROUTINE SIMULATES A VARIABLE BANDWIDTH y 4--P0LE 
CHEBYCHEV LOW PASS FILTER. (SUBROUTINE DIOFIL DOES THE 
ACTUAL FILTERING. A---TIME SERIES DATAy ND=»# OF DATA POINTS y 
EFFICIENT FILENAME 


•NAME-H-ILTER COf 


BYTE ANSyYEAyFNAME(I5) 

LOGICAL FIRST 

D I MENS I ON A ( ND y 2 > y CR ( 8 y 5 ) y C I ( 8 y 5 ) y AR ( ND ) y A I ( ND ) 
DATA YEA /'Y'/ 

WRITE (50) A 


SELECT FILTER BANDWIDTH 


WRITE (7y200) 

FORMAT ('OFILE DESCRIPTION V 

8/0TEST:L CHEBYCHEV y 4-POLE y 0.5 DB RIPPLE y FCO-"IV 

%' TEST2 CHEBYCHEV y 4-POLE y 0.5 DB RIPPLEy FC0--:>.5'/ 

S'OINPUT FILTER COEFFICIENT FILENAME' ) 


READ FILTER COEFFICIENTS FROM FILE 


READ ( :L0 y 208 ) FCO y FS y RIP y NP y F'STP y HSTP y N y ( ( CR ( J y I ) y J=:l y 8 ) y I“-.l y N > 


CALL GETSTR(SyFNAMEy 14y lERR) 

CALL ASSIGN (10 y FNAME y 14 y ' RDO ' ) 

'Si 

FORMAT (/3(28XyG13.A/) y28XyI3/y2(28XyG13.6/) y33XyI2y//y 
8(2Xy8(2XyG13.6) ) ) 

CALL CLOSE (10) 

DO 20 J™ly8 
DO 10 I-"lyN 
CI(JyI)::=CR(JyI) 

CONTINUE 

CONTINUE 


PERFORM FILTERING VIA SUBROUTINE DIGFIL. DATA IS PASSED 
THROUGH FILTER TWICE TO OBTAIN INITIAL CONDITIONS. 


FIRST™. TRUE. 

C A L L D I G !"• I L ( 2 y C R y N D y A R y A I’! ) 

C AI..L D I Gl" I L. ( 2 y C I y ND y A I y A I ) 
IF(FIRST .EO.. FALSE.) GO TO 40 
FIRST™. FALSE. 

REWIND 50 
READ (50) A 
GO TO 30 
REWIND 50 
RETURN 
END 


I 
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